Absence of aging in the remanent magnetization in Migdal-Kadanoff spin glasses 
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We study the non-equilibrium behavior of three-dimensional spin glasses in the Migdal-Kadanoff 
approximation, that is on a hierarchical lattice. In this approximation the model has an unique 
ground state and equilibrium properties correctly described by the droplet model. Extensive nu- 
merical simulations show that this model lacks aging in the remanent magnetization as well as 
a maximum in the magnetic viscosity in disagreement with experiments as well as with numerical 
studies of the Edwards- Anderson model. This result strongly limits the validity of the droplet model 
(at least in its simplest form) as a good model for real spin glasses. 
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Spin glasses are disordered magnets which for low im- 
purity concentrations above the Kondo regime display 
interesting non-equilibrium phenomena. In particular, a 
freezing of the dynamics appears at a temperature T c be- 
low which slow relaxation phenomena manifest through 
non-stationary effects in the zero-field-cooled magnetiza- 
tion. In this regime different non-equilibrium phenom- 
ena have been observed such as aging, remanence and 
several memory as well as chaotic effects pj . Despite the 
great activity devoted to understand the nature of the low 
temperature phase in three-dimensional spin glasses (nu- 
merical simulations, experiments and theory, see the re- 
view @]) still many questions regarding the ground state 
(e.g. its shape and its uniqueness) and the type of exci- 
tations remain unanswered. 

The mean-field picture for spin glasses || (i.e. the re- 
sults obtained for the Shcrrington-Kirkpatrick model), 
despite its great theoretical interest, it is not able to fur- 
nish a real space picture of the type of excitations present 
in spin glasses. To fill this gap, and based on domain 
wall scaling arguments (initially proposed by McMillan 
and Bray and Moore), Fisher and Huse proposed what 
has been termed as droplet model for spin glasses [H . In 
the droplet model there are two unique ground states re- 
lated by spin inversion symmetry. Thermal fluctuations 
activate droplets which are supposed to be compact do- 
mains of typical size L and fractal surface of dimension 
d s > d — 1. These excitations cost a free energy which 
grows like T(T)L where 8 is a zero-temperature expo- 
nent and T(T) is a temperature dependent stiffness con- 
stant. The idea that excitations in spin glasses are com- 
pact droplets is the simplest description that finds its 
most successful application in the study of phase transi- 
tions in ordered systems. Despite its inherent simplicity, 
the droplet model has a severe limitation, i.e. its main 
assumptions remain to be proven from a correct micro- 
scopic theory. If one of its key assumptions were wrong 



then the whole set of predictions coming out from the 
model should need to be revisited. 

Now, what is an appropriate microscopic model which 
describes the spin glass transition? The simplest pro- 
posal was putted forward by Edwards and Anderson al- 
most twenty years ago B, who introduced a random 
bond nearest-neighbor interaction model, the so called 
Edwards- Anderson model (EA model). It is widely be- 
lieved that the EA model is a real spin glass, i.e. it repro- 
duces the major part of results experimentally measured 
in the laboratory. So the question is whether the droplet 
model is the appropriate theory to describe the phe- 
nomenology already contained in the EA model. Despite 
of the large number of numerical works devoted to this 
question (see the reviews MM) there is still no an uni- 
versal agreement on it. 

Our work has been motivated by recent results by the 
Manchester group || who found that finite-size effects 
in the Migdal-Kadanoff approximation (MKA) of the 
three-dimensional EA model are mean-field like. While 
in the thermodynamical limit the MKA is known to be 
described by the droplet model with d s = d — 1 and 
6 ~ 0.26 pJlQl. So the Manchester group suggested that 
the droplet model could also explain the vast majority of 
numerical simulation results for the EA model obtained 
during the last decade (which on the other hand, have 
been taken by the Rome group as evidences against the 
droplet picture). This is an interesting observation whose 
physical meaning and consequences need to be better 
understood and was already anticipated quite long ago 
in a theoretical study of the one-dimensional spin-glass 
chain pd| . This controversy has been centered around 
the study of the spin-glass equilibrium properties. So, it 
is now time to check whether non-equilibrium behavior 
is well reproduced by the droplet model. This is of the 
outmost importance because experimental measurements 
in spin glasses in the low temperature regime are always 



taken in the out-of-equilibrium regime. 

In this paper we want to show that the droplet model 
lacks one of the key features of real spin glasses found 
in the laboratory, i.e. the aging in the zero-field-cooled 
magnetization. Consequently, the physics contained in 
the droplet model corresponds to a limited class of dis- 
ordered systems being far from what is observed in real 
spin glasses. 

The EA model in the presence of a field is defined by, 



T~t = — 2_^ JijViVj — h 2_^ &i 



(1) 



where the site indexes run on the nodes of a cubic lat- 
tice, (i,j) stands for nearest neighbor pairs, the spins take 
values <7j = ±1 and the couplings are extracted from a 
Gaussian distribution of zero average and unit variance. 
Following H we will consider the three dimensional EA 
model in the MKA which amounts to consider a hierar- 
chical lattice that is constructed iteratively by replacing 
each bond by eight bonds as indicated in fig. [y. Denoting 
by g the number of generations then the total number of 
bonds is 8 9 which corresponds to the number of sites for 
a cubic lattice with lattice size L = 2 9 . 




FIG. 1. Elementary step in the construction of the hierar- 
chical lattice, where the MK approximation is exact. 

The order parameter can be defined through the equi- 
librium autocorrelation function, 



= lim lim 

t— >0O V^cx 



Ei=i_gjM^iW) 



(2) 



where V — 8 9 — L 3 is the volume and the averages (..) 
are taken over dynamical histories starting from differ- 
ent equilibrium initial conditions at time 0. The pa- 
rameters Xi are weights which may consider the fact 
that a given site is connected with a different number 
of bonds depending on its generation level (i.e. depend- 
ing on which iteration in the recursive construction of the 
lattice that site was generated). Our results concentrate 
on the choice x% — 1, i.e. all sites are identically weighted. 
However also the results obtained with Xi = Ci, where Cj 
is the connectivity of site i (so all bonds are identically 
weighted) corroborate our conclusions |12| . 

We have concentrated our attention in the study of 
the relaxational dynamics in the low-temperature phase 
T <T C ~ 0.88 H below which q EA defined in Eq.(|) is 
different from zero. We have used Monte Carlo dynamics 
with Metropolis algorithm and random updating |13J . In 
our runs we follow the typical aging experiment schedul- 
ing, that is: at t = we quench the system from infinite 
temperature to a finite one T < T c without magnetic 



field, letting the system evolve for a time t w . At time t w 
we switch the magnetic field on. For subsequent times 
(t > t w ) the system continues to relax in a magnetic field 
h and then we measure the following two quantities: a) 
the autocorrelation function 



G (t, t w ) 



J2i=l X i a i( t ) CT i( t -» 



(3) 



and b) the zero-field-cooled susceptibility defined by 

XzFC(t,t w ) = hm , 4 

h^O h 



where 
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l^i=i Xi 



<Ji{t) - ai(t w )] 



Z-*ii=l X i 



(5) 



The limit in cq.(]J) is usually ignored, because we al- 
ways work in the linear response regime. All the data 
we present have been obtained with a magnetic field of 
intensity h — 0.1 and we have checked that the same sus- 
ceptibility is obtained by doubling the perturbing field. 

We have performed extensive numerical simulations for 
,g = 5 (L = 32) and g = 6 (L = 64) at two different tem- 
peratures (T = 0.7,0.5) and for many values of t w . We 
obtain the same results for both temperatures. Here we 
present only those for T = 0.7, while the complete set 
of data will be reported elsewhere. Note that the ratio 
T/T c ~ 0.8 used here is very similar to that used in many 
experiments |l4[ . 

In fig. we show the autocorrelation function for dif- 
ferent values of t w . One can observe the presence of aging 
characteristic of many glassy systems. 
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FIG. 2. The correlation functions for g = 6, T = 0.7 and 
different t w clearly show aging. In the inset the same curves 
as a function of the scaling variable log(i)/log(iu,). 

Following general assumptions, in the asymptotic 
regime t w — > oo, the correlation function decomposes in 
two terms each one governing a different time regime. 
In the quasi-equilibrium regime t — t w << t w the sys- 
tem is in some sort of local equilibrium and correlation 



functions are time translational invariant. In the aging 
regime t—t w » t w the system ages and correlation func- 
tions depend on both times through non trivial scaling 
relations. So in general one can write [|15fl, 



L^[Z 1 T W ) — L^ s t[t t w ) ~r L; a gi n g [T,T W ) 



(6) 



with linv^oo C s t(r) = q EA . In equilibrium the aging part 
vanishes and one recovers the previous result of eq. (g) . 

The very difference between the experimental data and 
the EA model on one hand and the MKA and the droplet 
model on the other is the large times scaling of dynami- 
cal functions. As can be seen in the inset of fig. g, in the 
MKA we find that the aging part of the autocorrelation 
function is well described, in the large times limit, by a 
function of the ratio \og(t) / \og(t w ) . On the other hand in 
experiments and in the EA model the scaling is far from 
the \og(t) / \og(t w ) and similar to t/t w . Because of the 
use of the scaling variable, in the inset of fig. g the data 
corresponding to the quasi-equilibrium regime collapse 
on the line \og(t) / \og(t w ) — 1. Thus we can estimate 
the value of q EA as the limit of the scaling function for 
log(t) / \og(t w ) — » 1+, i.e. q EA ~ 0.6 (a value compatible 
with §]). 
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FIG. 3. The zero-field-cooled susceptibility data do not 
show any aging. In the legend we report the values of g and 
tw In the inset we show the same quantity measured in the 
EA model, which resembles a lot the experimental data. 

Our most striking result is found for the zero-field- 
cooled susceptibility XZFc(t,t w ) shown in fig. g. In the 
MKA there is no dependence of the susceptibility on t w . 
We believe that such a result, which is characteristic of 
droplet models and kinetic growth, makes the droplet 
model, at least in its simplest form, inadeguate for the de- 
scription of the EA one. Aging in both zero-field-cooled 
and field-cooled magnetization is so commonly found in 
experiments on spin glasses that it is not clear to us how 
this result can be explained by the standard droplet the- 
ory. Note also that the peak in the magnetic viscos- 
ity S(t,t w ) — dxzFc{t,t w )/dlog(t-t w ) (experimentally 
very well observed Jig) is completely absent in the MKA. 



We should remind that aging in XZFc(t, t w ), with a peak 
in the S(t,t w ), is naturally found in the EA model (see 
inset of fig. ||) as well as in mean-field models. Then it 
remains to be explained why these aging effects are nat- 
urally and easily observed in the EA model and not in 
the MKA. 

Finally we consider the analysis of the fluctuation- 
dissipation ratio useful to compare the results obtained in 
the MKA with those obtained in the EA and coarsening 
models |17| ] . In the quasi-equilibrium regime (t — t w << 
t w ) the system is in local equilibrium. Consequently both 
correlation and susceptibility are time-traslational invari- 
ant and the fluctuation-dissipation theorem (FDT) is sat- 
isfied, 



T X ZFc(t -t w ) = l-C(t- t w ) 



(7) 



In the aging regime (t — t w » t w ) the system ages and 
FDT is violated. Then it is useful to define the so called 
fluctuation-dissipation ratio |18| 



X(t,t w ) 



TR(t, t w ) 

dC(t,t w ) 



(8) 



which, in the asymptotic long-times limit t, t w — ^ oo, 
may be uniquely expressed as function of the correlation 
C{t,t w ) yielding 

T X ZFc(t,t w )= f X(C)dC . (9) 

JC(t,t w ) 

Moreover the X(C) can be related to equilibrium quan- 
tities fl9[j . The previous expression reduces to Eq.(J7]) 
in the quasi-equilibrium regime where X = 1. A plot 
of TxzFc{t,t w ) as a function of C(t,t w ) is expected to 
show two different behaviors. For q EA < C < 1 we have 
X = 1 and so the curve Txzfc versus C has slope — 1. 
For C < q EA the X may be a non- vanishing function of C 
andwehaveT XZ Fc(t,t w ) = {1 - q EA ) + jf {t A tw) X (C)dC . 
In coarsening models, X = for C <q EA and so the func- 
tion xzFciC) is flat for C < q EA . In fig. || we show 
the xzfc as a function of C for different values of g 
and t w , which show that the behavior rapidly converges 
to that of coarsening models and strongly differs from 
that observed in finite-dimensional EA spin glasses [|L7| . 
The horizontal line in fig. is the infinite time limit of 
the susceptibility, extrapolated from the data of fig. || 
and from those for the field-cooled magnetization (not 
shown). It is an upper bound for the plotted curves, 
thanks to the positiveness of X ratio. From fig. H we can 
also get an estimate for the q EA order parameter, defined 
as the abscissa value where the curves leave the FDT line 
(Txzfc = 1 — C). Very reasonably this point is converg- 
ing, in the large times limit, near to the intersection of 
the two lines, giving q EA ~ 0.6 (as already found from 
the data of fig. ||). 







discussion on finite-dimensional spin glasses. 
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FIG. 4. The absence of response in the aging regime is sug- 
gested by the rapid convergence of the xzfc versus C curves 
to the plotted lines. Note that the horizontal line is an upper 
bound for the data. 

Fig. [| adds more evidence on the fact that spin glasses 
in the MKA do not capture all the key features of finite- 
dimensional spin glasses as we know them from the three- 
dimensional EA model. 

To summarize, we have shown that in the MKA spin 
glasses do not show aging in the integrated response func- 
tion. This aging is experimentally observed in real spin 
glasses through zero-field-cooled and field-cooled mea- 
surements being one of the key features which distin- 
guishes spin glasses from other disordered systems. The 
study of the fluctuation-dissipation ratio suggests that re- 
laxation in this model is driven by coarsening, like in con- 
ventional ferromagnets. One could argue that these re- 
sults for the MKA are not extensible to the droplet model 
because, in the general case, the inequality d s > d — 1 
could restore aging. Despite this possibility our results 
unambiguously show that the MK model is not a good 
model for realistic spin glasses. A new class of excitations 
or droplets must be present in spin glasses. The droplet 
model in its simplest version does not capture the physics 
behind real spin glasses. 

One possible generalization of the droplet model 
(which would be no longer simply droplet) is to consider 
two kind of basic excitations in a spin glass: on small 
length scale the usual droplets and, in addition, system- 
size scale collective rearrangements p(| • The second kind 
of excitation are, at present, ignored in the droplet model 
(they are exponentially rare), but they could be respon- 
sible for the many mean-field like features observed in 
finite-dimensional spin glasses. In terms of a very sim- 
plified energy landscape the two excitations would corre- 
spond, respectively, to the local movements of the system 
in a single "valley" and to the jumps from one valley to 
an other one. In our opinion a new theory comprehen- 
sive of the small scale droplets and the system-size scale 
excitations (with a clear real space picture) would be 
welcome and could hopefully terminate the longstanding 
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